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ABSTRACT 

We present the results of a combined X-ray and Ha study of 10 galaxy groups and 17 galaxy 
clusters using the Chandra X-ray Observatory and the Maryland Magellan Tunable Filter. We find 
no difference in the morphology or detection frequency of Ha filaments in groups versus clusters, over 
the mass range 10^^ < M500 < lO^^M©. The detection frequency of Ha emission is shown to be 
only weakly dependent on the total mass of the system, at the 52% confidence level. In contrast, we 
find that the presence of Ha filaments is strongly correlated with both the global (89% confidence 
level) and core (84%) ICM entropy, as well as the X-ray cooling rate (72%). The Ha filaments are 
therefore an excellent proxy for the cooling ICM. The Ha filaments are more strongly correlated with 
the cooling properties of the ICM than with the radio properties of the BCG; this further supports the 
scenario where these filaments are directly associated with a thermally-unstable, rapidly cooling ICM, 
rather than radio bubbles. The ICM cooling efficiency, defined as the X-ray cooling rate per unit gas 
mass, is shown to correlate with the total system mass, indicating that groups are more eflacient at 
cooling than clusters. This result implies that, in systems with cool cores, AGN feedback scales with 
the total mass of the system, in agreement with earlier suggestions. 

Subject headings: galaxies: cooling fiows - galaxies: clusters - galaxies: groups - galaxies: elliptical 
and lenticular, cD - galaxies: active - ISM: jets and outfiows 



1. INTRODUCTION 

The high X-ray surface brightness in the cores of some 
galaxy clusters suggests that, in the absence of any feed- 
back, radiative cooling should become a runaway pro- 
cess, leading to "cooling fiows" of 10-100 M© yr~^ (see 
review by Fabian 1994). A distinguishing characteris- 
tic of cool core clusters is the presence of warm, optical 
line-emitting gas (Hu et al. 1985, Heckman et al. 1989, 
Crawford et al. 1999, Jaffe et al. 2005, Hatch et al. 2007). 
Typically observed at Ha, this warm gas generally ex- 
hibits either a nuclear or filamentary morphology and 
has been shown to correlate with several properties of 
the cool core, such as cooling rate and cluster entropy 
(McDonald et al. 2010; hereafter M-hlO). Despite the 
near ubiquity of optical line-emitting nebulae in cool core 
clusters, their origins remain unclear. In M+10 we con- 
sidered several processes for producing the observed Ha 
morphology, and found that the observations were con- 
sistent with two formation scenarios: 1) the warm gas 
cooled out of the X-ray, and 2) the warm gas was en- 
trained behind buoyant radio bubbles arising from AGN 
feedback (e.g. Churazov et al. 2001, Reynolds et al. 2005, 
Vernaleo et al. 2007, Revas et al. 2008). The high Ha 
surface brightnesses measured in these systems indicate 
that additional ionization sources are present, allowing 
the ionized gas to recombine more than once. Crawford 
et al. (2005) describe potential sources of ionization, the 
most promising of which are: 1) cosmic ray ionization, 
2) heat conduction from the ICM to the colder filaments, 
and 3) photoionization by hot young stars. In McDonald 
et al. (2009) and M+10 we provide evidence for both star 
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formation and conduction, respectively, but were unable 
to rule in favor of a dominant ionization mechanism. 

While much effort has gone into quantifying and ex- 
plaining the presence of warm filaments in clusters, very 
little work has been done on the environmental depen- 
dence of this phenomenon. While massive clusters of 
galaxies are dominated by gravitational and thermal 
processes, their low-mass counterparts, galaxy groups, 
may be dominated by non-gravitational processes such as 
winds, cooling fiows or AGN feedback (Sun et al. 2009a). 
Thus, by considering low-mass systems we can hope 
to find out whether the observed phenomena are more 
closely related to global, gravitational processes or more 
local, baryon physics. In M+10, we showed that the 
properties of the Ha filaments were strongly correlated 
with the properties of the cluster core, such as the cooling 
rate, cooling radius, entropy and temperature. However, 
most of these trends only cover a small range in X-ray 
properties, since M+10 only considered clusters. 

In an effort to properly evaluate the dependence of 
the Ha filaments on environment, we have conducted a 
survey of 10 galaxy groups using the Maryland Magellan 
tunable filter (hereafter MMTF; Veilleux et al. 2010). We 
combine this sample of 10 groups with a subset of 17 clus- 
ters drawn from M+10. In §2, we describe this sample in 
more detail, along with the data acquisition and analysis. 
In §3, we present the results of this study, re-examining 
several trends found by M-hlO, and also exploring new 
issues. In §4 we discuss the implications of these results 
and postulate on the origin and power source of these 
warm filaments. Finally, in §5 we conclude and discuss 
prospects for future studies. 

Throughout this paper, we assume the following cos- 
mological values: Hq = 73 km s~-^ Mpc~-^, ^matter = 

0.27, ^vacuum ~ 0.73. 
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2. DATA COLLECTION AND ANALYSIS 
2.1. Sample Selection 

The goal of this study is to extend our investigation of 
Ha filaments (M+10) from the cluster to group environ- 
ment. Thus, we begin with a sample drawn from M+10, 
of 17 cool core clusters. These clusters were originally 
chosen from the larger sample of White et al. (1997). 
This large sample, much like any X-ray selected sam- 
ple, is inherently biased towards cool cores due to their 
higher surface brightness. The full sample of 207 clusters 
(White et al. 1997) was reduced based on the following 
criteria: 1) visible with the Magellan telescope {S < 35°); 
2) appropriate redshift to be imaged at Ha with MMTF 
(0.0 < z < 0.09). From this reduced sample, we choose 
17 clusters with high-quality Chandra imaging, covering 
three orders of magnitude in classical cooling rate. These 
clusters were chosen to yield a relatively flat distribution 
of cooling rates, ensuring that we covered the full gamut 
of cool core and non cool core clusters. This selection 
technique will, of course, introduce a significant bias, the 
effects of which will be discussed later in §4. The 17 se- 
lected clusters all have deep Ha imaging from MMTF, 
as presented in M+10. 

From the sample of 43 galaxy groups defined by Sun 
et al. (2009a), we chose at random 10 groups which obey 
the same criteria as above, with average temperatures 
ranging from 0.7-3 keV, to add to this sample. Unlike the 
White et al. (1997) sample, this collection of 43 galaxy 
groups is not based solely on X-ray detections of clusters. 
Sun et al. (2009a) include AGN-selected and optically- 
selected groups, removing much of the bias towards cool 
cores that affects X-ray selected catalogs. Thus, the 10 
groups in our sample should be relatively unbiased in 
mass, temperature and cooling rate. 

The full sample is listed in Table |5| We discuss below 
the processing of the new MMTF Ha observations and 
the archival Chandra X-ray data. 

2.2. Ha: MMTF 

MMTF has a very narrow bandpass (~5-12A) which 
can be tuned to any wavelength over ~5000-9200A 
(Veilleux et al. 2010). Coupled with the exquisite image 
quality at Magellan and the wide field of the Inamori- 
Magellan Areal Camera & Spectrograph (IMACS), this 
instrument is ideal for detecting emission-line filaments 
in galaxy groups and clusters. In M+10 we presented 
deep MMTF Ha observations of 23 clusters (17 of which 
have high-quality CXO data). As a follow-up to that 
project, we have observed an additional 10 groups during 
2009-10 at both Ha (A=6563A) and continuum (± 60A), 
for a total of 40 minutes each with the largest available 
bandpass 12 A). The typical image quality for these 
exposures was 0.6 ±0.2'^ These data are a significant im- 
provement on previous narrow-band imaging of cluster 
cores due to the very narrow (~ 10 A) bandwidth of the 
MMTF. This allows us to isolate and measure the flux of 
the Ha line without making any assumptions about the 
[N n]/Ha ratio. 

These new data were fully reduced in exactly the same 
way as the M+10 data, using the MMTF data reduction 
pipeline^, which performs bias subtraction, flat fielding, 
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sky line removal, cosmic ray removal, astrometric cal- 
ibration and stacking of multiple exposures (following 
Veilleux et al. 2010; see also Jones, Bland-Hawthorn, 
Shopbell 2002). The continuum image was then PSF 
and intensity matched to the narrow-band images to al- 
low for careful continuum subtraction. The stacked im- 
ages were calibrated using spectrophotometric standards 
from Oke (1990) and Hamuy et al. (1992, 1994) The error 
associated with our absolute photometric calibrations is 
~15%, which is typical for tunable filters and spectro- 
graphs. Finally, the data were corrected for Galactic 
extinction, following Cardelli et al. (1989) using redden- 
ing estimates from Schlegel et al. (1998). We do not 
attempt to correct for intrinsic extinction since the dust 
content of the optical filaments is not well known. All 
of these procedures are described in detail in Veilleux 
et al. (2010). 

For systems with complicated morphologies. Ha fluxes 
were measured by creating (by eye) a region which gen- 
erously traced the Ha emission and calculating the total 
signal within this region. For more symmetric morpholo- 
gies, a circular aperture centered on the emission peak 
was used, with the radius chosen to contain all of the 
obvious emission. We show in M+10 that this technique 
yields fluxes that agree relatively well with those in the 
literature. 

2.3. X-Ray: Chandra 

Archival data from the Chandra X-ray Observatory 
were retrieved for all 10 of our sample groups. In order 
to ensure a homogeneous treatment of the X-ray data, 
the clusters data were re-reduced alongside the groups 
data. These data were reprocessed with CIAO (version 
4.1.2) and CALDB (version 4.1.1) using the latest time- 
dependent gain adjustments and maps to create new level 
2 event files. Due to the large angular extent of some 
of the groups and clusters in our sample, we were re- 
quired to construct blank-sky background event files, us- 
ing the ACIS blank-sky background database, to prop- 
erly account for background flux. The new level 2 event 
files were cleaned for flares, using the Ic.clean routine, 
by examining the light-curve and removing any spurious 
bursts in intensity. These data cleaning and calibration 
procedures are all outlined in detail in the CIAO science 
threads^ . 

In order to separate the filaments or other substruc- 
tures from the X-ray halo, we applied an unsharp mask 
to each image, subtracting a 10'^ Gaussian smoothed im- 
age from a 1.5^' Gaussian smoothed image. The resulting 
image highlights any fine structure in the X-ray morphol- 
ogy (see Figure [T]). 

For each object, background-subtracted spectra were 
extracted using dmextract. Updated response files were 
created using mkacisrmf and mkwarf^ following the 
CIAO science threads. Counts were grouped into bins 
with 20 counts per bin, over the range 0.3 to 11.0 keV. 
Spectra were extracted in a variety of regions to better 
understand the relationship between the ICM and the 
Ha emission. These regions include: i) circular annuli, 
with spacing chosen so that Tout/^in = 1.25-1.6 (follow- 
ing Sun et al. 2009a, McDonald et al. 2010), ii) Coinci- 
dent with the Ha-emitting filaments, and iii) the regions 

^ http:/ /cxc. harvard.edu/ciao/threads/ 
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surrounding, but not overlapping with, the Ha-emitting 
filaments. 

In order to derive physical quantities from the X-ray 
spectra, we used the XSPEC spectral fitting package (Ar- 
naud 1996). For the circular annuli, we first deproject 
the data in each radial bin using the direct spectral de- 
projection method (DSDEPROJ; Russel et al. 2008). We 
then model the spectra with a combination of photoelec- 
tric absorption (PHABS) and thermal brehmsstrahlung 
emission from a hot diffuse gas (MEKAL). This combi- 
nation of models has 3 free parameters which describe 
the state of the ICM: Tx, ng, and Z. In order to mea- 
sure the strength of the ICM cooling, we also allow for 
an additional component which represents gas cooling 
over a range of temperatures (MKCFLOW), which has 
an additional free parameter, dM/dt. 

For a more complete description of the spectral extrac- 
tion and modeling techniques, the reader is directed to 
McDonald et al. 2010. 

3. RESULTS 
3.1. Warm Ionized Filaments 

The primary strengths of MMTF are its large field 
of view and excellent angular resolution. These make 
it ideally suited to search for thin, extended filaments 
in the cores of clusters, as we showed in M+10. Our 
previous study of galaxy clusters revealed a wide vari- 
ety of filament morphologies, leading us to develop a 
classification scheme to help with the analysis of such 
a diverse sample. We found that 35% of clusters ex- 
hibited thin, extended Ha filaments (type I), 30% had 
nuclear or only marginally extended Ha emission (type 
II), and 35% had no detectable Ha emission at all. We 
adopt this same nomenclature for our new sample of 10 
galaxy groups, finding very similar ratios between the 
three types (30% : 30% : 40%). The similar frequency 
with which we detect Ha filaments in groups compared 
to clusters suggests a similar formation mechanism. 

In Figure [T] we show the new data for 10 galaxy groups. 
For all 10 groups in this sample we show the MMTF 
red continuum, smoothed CXO, unsharp masked CXO, 
and MMTF Ha images. Just as we saw for clusters, the 
presence of complex morphology in the X-ray appears 
to correlate with the presence of Ha filaments. We find 
that all three clusters with strongly asymmetric residuals 
in the unsharp masked X-ray images have extended Ha 
filaments. We will investigate this relationship between 
the X-ray and Ha data further in the following sections. 

Of particular interest in Figure [l] is the morphology 
of the Ha emission for Abeh 1991 and NGC 4325 (Ser- 
sic 159-03 was discussed in M+10). Abeh 1991 has 
an arrow-shaped morphology, extending north from the 
BCG nucleus. This morphology is reminiscient of a bow- 
shock. We find that the arrow-head corresponds to a 
bright blob of soft X-rays, while the peak of Ha emis- 
sion at the base of the arrow has no X-ray counterpart. 
The second interesting group, NGC 4325, exhibits sev- 
eral distinct, radial filaments. This complex morphology 
strongly resembles the core of the Perseus cluster. Un- 
like Perseus, the central galaxy in NGC 4325 is almost 
completely isolated and the total X-ray mass is an order 
of magnitude lower. This suggests that the global prop- 
erties of the group/cluster plays only a minor role in the 



presence and morphology of Ha emission. 

3.2. Global and Core X-Ray Properties 

In order to compare properties of groups and clusters, 
we adopt a scale radius at which we measure quanti- 
ties such as temperature, density, entropy and enclosed 
mass. The radius we use, r25oo, is the radius at which the 
average enclosed density is 2500 times the critical den- 
sity of the Universe. This radius is considerably smaller 
than the more commonly used value of rsoo, however we 
prefer using r25oo due to our ability to directly measure 
quantities at r25oo (which we will refer to as "global" pa- 
rameters), rather than infer them indirectly from other 
measurements. The typical value of rsoo for the groups 
in our sample is ~ 600 kpc (Sun et al. 2009a), which cor- 
responds to a diameter of ~ 20' at a typical redshift of 
0.05 - larger than the typical CXO field of view. Thus, 
we choose a smaller radius at which we can measure var- 
ious properties for all of the groups and clusters in our 
sample. 

In Figure |2] we plot various X-ray scaling relations for 
our full sample. For systems in hydrostatic equilibrium, 
we expect a direct correlation between the total enclosed 
mass and the gas temperature for a given radius, which 
we observe in general. However, we also find several sys- 
tems which deviate from this relation. This is because 
clusters are often not fully in hydrostatic equilibrium due 
to processes such as AGN feedback or mergers. Never- 
theless, these scaling relations show the broad range in 
parameters covered by the groups and clusters in this 
sample. 

At a glance. Figure |2] offers no new insight into the 
presence of Ha filaments. We detect Ha in systems at 
all temperature, mass and entropy, with no obvious bias. 
To further quantify this, we consider both the fraction of 
systems with non-zero Ha fiux and those with extended 
Ha filaments as a function of global X-ray properties 
in Figure 13] We have chosen the binning such that the 
number oisystems in each bin are roughly equivalent. In 
each panel, we show the F-test likelihood that the his- 
togram is fiat, compared to monotonically increasing or 
decreasing. This statistic confirms that, at the 87% con- 
fidence level, the presence of Ha emission is uncorrelated 
with the global temperature. Additionally, we find only 
a weak correlation between the presence of Ha emission 
and system mass (52% confidence level). 

We note that, in higher mass/temperature systems, the 
fraction of systems with Ha emission is likely an upper 
limit, due to our bias towards cool cores in these sys- 
tems. Thus, while there is only a weak correlation be- 
tween the presence of Ha emission and the global mass, 
the trend would likely be strengthened by the inclusion of 
high mass, non cool-core systems which we would not ex- 
pect to emit at Ha based on our previous work (M+10). 

The bottom two panels of Figure show the distribu- 
tion of the entropy and the gas mass fraction for systems 
with non-zero Ha fiux or with extended Ha filaments. 
We find a correlation between entropy and the presence 
of Ha filaments (89% confidence level). We find that 
12 of the 15 systems with K < 900 keV cm^ show fil- 
aments, while only 5 of the 11 systems with K > 900 
keV cm^ show filaments. This is consistent with findings 
by previous studies (Cavagnolo et al. 2008, McDonald 
et al. 2010) which show that low-entropy systems tend 
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to have more star formation and optical line emission. 
Finally, we observe very few systems with Ha filaments 
and very low X-ray gas mass fractions, suggesting that 
the two gas phases are related (51% confidence level). 
This seems to favor the hypothesis that the warm ion- 
ized gas is intimately linked to the hot gas - a certain 
threshold amount of hot gas is needed in order for Ha 
filaments to exist. In summary. Figure [3] suggests that 
the presence of Ha emission is only weakly dependent on 
the global mass (and thus, temperature) of the system. 
However, the state of the ICM (i.e. entropy, fractional 
gas mass) may help dictate whether these filaments can 
exist. 

In M+10, we provided evidence that the observed Ha 
filaments were intimately linked to the X-ray cooling. 
Thus, we re-examine the properties of the group /cluster 
cores, where the gas is cooling on a short timescale. 
Figure [4] provides similar histograms to those described 
above, but now considering only measurements taken in 
the group/cluster core. We calculate the frequency of 
Ha emission as a function of the 1.4 GHz radio power, 
which traces AGN feedback, the X-ray cooling rate, the 
average temperature within the central 100 kpc, and the 
entropy at a radius of 50 kpc. We find a weak correla- 
tion between the 1.4 GHz luminosity and the presence 
of Ha filaments (56% confidence level). In contrast, the 
correlation with the X-ray cooling rate is significantly 
stronger (72%), as we also found in M+10. This sug- 
gests that the ionized filaments are linked more closely 
to the cooling ICM than to the radio-loud AGN. In the 
lower panels, we see a much stronger dependence on tem- 
perature and entropy when we move to smaller radius. 
There appear to be very few systems with Ha filaments 
and high temperature/entropy in the inner regions. This 
confirms our findings from M+10 that the presence of 
Ha filaments correlates with the temperature (55% con- 
fidence level) and entropy (84% confidence level) inside 
of the cooling radius. The results of Figure [4] are not 
strongly affected by our selection biases. The missing 
clusters, which should be high-mass with no cool core, 
should have high entropy and temperature in the inner 
100 kpc. Since we expect these non-cooling systems to 
lack Ha emission, the trends would remain the same. 

The combination of Figures [3]-[4] suggests that the pres- 
ence of Ha filaments is a function of the state of the ICM 
inside of ~ rcooi and is largely independent of the global 
properties much beyond this. 

3.3. X-Ray - Ha Correlations 

In M+10 we provided several pieces of evidence which 
link the observed Ha filaments to the X-ray cooling. 
With the addition of 10 additional low-mass systems, it 
is relevant to return to these results and ensure that they 
are still significant. 

In Figure [5] we plot both the total Ha luminosity and 
the Ha luminosity in filaments as a function of various 
X-ray properties of the cool core. In M+10 we found that 
systems with warm cores (/cTioo > 4.5 keV) do not emit 
at Ha, while systems with cool cores can have Ha emis- 
sion, but not always. This trend is preserved by includ- 
ing cooler systems (blue points). Additionally, we saw 
correlations between the Ha luminosity in filaments and 
the X-ray cooling rate and core entropy values. These 
trends are strengthened by the addition of low-mass sys- 



tems, which exactly follow the distribution of high-mass 
systems. The strong correlation between the Ha lumi- 
nosity and the X-ray cooling rate suggests a direct link 
between the cooling ICM and the warm gas. 

In addition to a correlation between the mass of gas 
cooling below the hot phase and the mass of warm gas, 
we also found, in M+10, a correlation between the ex- 
tent of the warm gas and the cluster cooling radius. The 
addition of 10 additional groups to this result (Figure 
[6| further strengthens our claim that the size of the Ha 
maments do not exceed the cluster cooling radius. Since 
the location of the cooling radius is dependent on an 
arbitrarily-assigned cooling time, we have tested cooling 
times of both 3Gyr and 5Gyr. While M+10 showed that 
a cooling radius based on tcooi = 5Gyr matched the data 
well, we show in Figure 6] that tcooi = 3Gyr does a sig- 
nificantly better job of defining the maximum radius of 
Ha emission. The fact that we do not see Ha emission 
beyond the radius at which the hot ICM is cooling in less 
than 3Gyr suggests a natural timescale for the formation 
of these filaments. 

In §2.2, we discussed the procedure of extracting X- 
ray spectra coincident to and surrounding the observed 
Ha filaments. Since we are unable to determine the 
3-dimensional shape of these filaments, we attempt to 
model them using two different geometries: 1) thick 
slabs of gas extending into the sky for the full length of 
the cluster, modeled with a single-temperature plasma, 
and 2) thin, cylindrical filaments that are modeled with 
a two-temperature plasma to account for the back- 
ground/foreground gas seen in projection. In Figure [7 
we show the results of this exercise for the groups anc. 
clusters in our sample which exhibit extended filaments. 
As we found in M+10, the X-ray gas coincident with the 
Ha filaments (labeled "in") appears to be cooling faster 
than the surrounding ICM (labeled "out"). Assuming 
thin-filament geometry, the X-ray temperature of the fil- 
aments is ~ 20-50% that of the off-filament gas at the 
same radius for most of the systems. Additionally, both 
the entropy and cooling time of the filaments is roughly 
an order of magnitude smaller for gas in the filaments 
than off-filament gas at similar radii. Even if we assume 
that the "filaments" are sheets seen in projection, we in- 
fer a cooling time half as long for in- filament gas. We 
find no significant difference in the filament properties 
between groups and clusters, suggesting that they share 
the same formation mechanism. 

Figure [7| offers unique evidence for a link between the 
hot and warm gas phases. The fact that the hot gas 
is cooling an order of magnitude faster in the filaments 
provides a direct connection between the X-ray cooling 
flow and the Ha filaments. Coupled with the correlation 
between the cooling rate and Ha flux contained in fila- 
ments, the evidence for a link between the Ha filaments 
and the X-ray cooling flow is further strengthened. The 
fact that we see no obvious differences between the prop- 
erties of the filaments in groups and clusters suggests 
that the formation mechanisms for these filaments are 
the same. We discussed formation scenarios in depth in 
M+10, but the addition of the lower mass systems now 
provide new insights in the discussion (next section). 

4. DISCUSSION 
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4.1. Groups vs Clusters: Differences and Similarities 

Despite the smooth transition in mass and temperature 
from groups to clusters, the former are not simply scaled 
down versions of the latter. We summarize the most 
relevant differences between the two below. 

• Feedback processes such as radio jets, starburst- 
driven winds, and merging with other bound sys- 
tems will dominate in the group environment, while 
they play second fiddle to gravitational processes in 
the cluster environment. At M < lO^^M© the total 
radio heating energy of the BCG AGN (assuming 
e = 0.1; Sun et al. 2009b) becomes larger than the 
total potential or thermal energy of the ICM. 

• While clusters form at late cosmological times {z < 
1), groups have been forming over almost the entire 
age of the universe {z < 10; Fakhouri et al. 2010). 
Since age roughly correlates with several X-ray ob- 
servables (entropy, gas fraction) one would expect 
groups to have a much broader range of properties 
due to their substantial spread in ages. 

• In galaxy clusters, the ratio of the total mass in 
gas to dark matter is roughly 1:10, with a negligi- 
ble contribution from light in stars using standard 
mass- to-light ratios. However, in galaxy groups, 
the stellar mass fraction is an order of magnitude 
higher, leading to significantly lower mass to light 
ratios (Giodini et al. 2009). Additionally, the gas 
fraction decreases as a function of total mass, such 
that groups can have stellar to dark matter ratios 
of roughly 1:10, with a negligible contribution from 
the ICM. Clearly this limits the influence that the 
ICM can have on the rest of the system. 

Remarkably, although there are such vast differences be- 
tween clusters and groups, we find very little dependence 
of Ha emission on the global mass or temperature of the 
group/cluster. The fact that we see similar cooling rates 
and optical line flux in systems with orders of magni- 
tude less gas implies that groups are more efficient in 
their cooling. In order to pursue this notion, we consider 
the cooling properties of the group/cluster as a function 
of the system mass (M25oo)- In Figure [s] we first show 
the X-ray mass deposition rate, dM/dt, as a function of 
the group/cluster mass, M25oo- In order to improve the 
mass coverage, we have added 6 additional groups from 
the Sun et al. (2009a) sample - these systems do not have 
Ha imaging. We observe a steady increase of the cooling 
rate as a function of total mass, which is not surprising 
since there is more gas available for cooling. In order 
to remove this dependence on system mass, we next con- 
sider the amount of gas cooling out of the X-ray, per solar 
mass of material in the group/cluster (upper right panel 
of Figure g. We see no observable trend between these 
two quantities. However, we do detect what looks like an 
upper limit to the cooling rate per unit mass, such that 
the maximum cooling is ~ 0.02% of the cluster mass per 
Gyr. 

If, instead, we consider the ICM cooling efficiency^ de- 
fined as the amount of gas cooling out of the X-ray per so- 
lar mass of gas in the group/cluster, we see a correlation 
(Pearson R=0.71). We can understand this correlation 



by considering the following: (1) The presence of cool 
cores is only weakly dependent on mass (Figure [s]); (2) 
The total amount of gas is strongly dependent on mass. 
Combining these two points naturally leads to a correla- 
tion between the ICM cooling efficiency and the system 
mass. Similarly, we can show the same trend using the 
Ha-derived star formation rate (directly proportional to 
Lffa)^ instead of the X-ray cooling rate, since these two 
quantities are strongly correlated. As we discussed in 
§2, this sample is biased against massive, non-cool core 
clusters. We point out that the correction of this bias 
would introduce additional high-mass, low effective cool- 
ing rate points to Figure [Sj As with Figure [Sj this cor- 
rection would act to strengthen the observed correlation. 
This correlation suggests that groups are more efficient 
at cooling than clusters, despite the fact that feedback 
processes dominate in this region. 

To better understand this trend, we can approximate 
the expected relationship of the ICM cooling efficiency 
with total mass: 
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Equation 1 makes the assumption that the local ICM 
density is roughly independent of the temperature, which 
we have confirmed for the groups/clusters in our sam- 
ple. If we couple equations 1 and 2 with the relation be- 
tween M and T found by Sun et al. (2009a; M oc T^ '^), 
we find the dependence of the ICM cooling efficiency on 
mass matches well with the observations (see Figure [8|. 
While this derivation is overly simplistic, the confirma- 
tion with observation is reassuring. A major issue with 
such a simple derivation is that it ignores the role of feed- 
back. In low-mass systems, one would expect the rela- 
tive energy contribution from the AGN to increase sub- 
stantially, yielding less efficient cooling. The fact that 
we observe the opposite suggests two possible scenarios 
for cool core clusters: i) AGN feedback plays a negligi- 
ble role in countering ICM cooling, or ii) AGN feedback 
scales with environment. The former explanation is in 
opposition to a rich literature on the subject, both the- 
oretical (e.g., Brighenti et al. 2003; Croton et al. 2006; 
Cattaneo et al. 2007; McCarthy et al. 2010) and obser- 
vational (e.g., Birzan et al. 2004; Donahue et al. 2006; 
Rafferty et al. 2008; Sun 2009b), which offers signifi- 
cant evidence for AGN feedback playing a dominant role 
in the heating of galaxy group/cluster cores. The sec- 
ond explanation, that AGN feedback scales with envi- 
ronment, is more plausible. This scenario was recently 
investigated by Sun (2009b), who found that, in cool core 
groups and clusters, the radio luminosity of the central 
AGN scales with the cooling luminosity. Sun found that 
cool cores in groups have consistently lower 1.4 GHz lu- 
minosity AGN than in clusters and suggested that the 
lower-pressure intragroup medium may not be able to 
sustain the radio-loud AGN seen in the cores of massive 
clusters. We elaborate on this idea in Figure [9| 
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Motivated by Figure 1 from Sun (2009b), Figure [9] 
shows the 1.4 GHz luminosity {Li^4ghz) as a function 
of the soft X-ray luminosity (I/o.5-2feey) inside of the 
cooling radius (Lo.5-2/cey)- We show in light grey points 
the distribution of groups and clusters from Sun (2009b), 
while the red ellipses show the loci of these distribu- 
tions. Sun identified two classes of system: coronae, 
which show no correlation between the cool X-ray and 
radio luminosities, and luminous cool cores, which have 
cool X-ray luminosity correlated with radio luminosity. 
We show in Figure [9] our proposed, simplified evolution 
of these systems. Luminous cool cores are experiencing 
feedback-regulated cooling, such that both the AGN and 
cool core are able to grow slowly over timescales on the 
order of Gyr, based on measured cooling times in the 
central lOOkpc. However, once the radio luminosity of 
the AGN reaches a certain threshold, it has enough en- 
ergy to equal the PdV work required to evacuate the cool 
gas from the central lOkpc - this threshold radio lumi- 
nosity is depicted by a dotted horizontal fine for three 
different mass regimes. Thus, systems experiencing suf- 
ficiently strong radio-mode feedback have disrupted cool 
cores where the radio and X-ray luminosities are decou- 
pled, yielding a corona-class object. Given the relatively 
few objects occupying the space between corona and cool 
core systems, the disruption of the cool cores must hap- 
pen relatively quickly (<cGyr). With no cooling fiow to 
fuel the central AGN, the radio luminosity will eventually 
drop until cooling is allowed again. Assuming the cool 
core is fully disrupted, it should take a few Gyr for gas 
which was originally at the cooling radius to form a new 
cool core. At this point, the central black hole becomes 
active once more and the feedback loop is re-established. 

In Figure |9] there are three systems with radio lumi- 
nosities sufficiently higher than the threshold for their 
mass - Hydra A, Abell 2052, and Abeh 4059. These 
systems are all known to be experiencing higher-than- 
normal radio-mode feedback and have severely disrupted 
cool cores (e.g. Wise et al. 2007; Reynolds et al. 2008; 
Blanton et al. 2009). These systems all lie slightly off 
of the locus of points for luminous cool cores, suggesting 
that they may be in the midst of evolving from cool core 
to corona class via a disruption of the cool core. 

In summary. Figures. [8]-[9] tell an interesting story 
about the cool cores of groups and clusters. The evidence 
suggests that groups are more efficiently converting the 
hot ICM into cool gas, as we detect Ha filaments and 
cool cores with similar luminosity in systems with orders 
of magnitude difference in gas mass. It is unsurprising 
that groups are cooling more efficiently, since the major- 
ity of the intragroup gas is already at low temperature. 
However, the fact that this cooling is allowed to proceed 
as expected suggests that AGN feedback, which we as- 
sume regulates the cooling process, scales with system 
mass as well. We confirm the findings by Sun (2009b) 
that this is indeed the case, with cool cores in groups 
harboring radio-quiet AGN, while cool cores in massive 
clusters tend to have radio-loud AGN. This scaling of 
feedback strength with environment is a result of cool 
cores in low mass systems being more easily disrupted 
than those in high-mass systems. 

4.2. The Origin of Ha Filaments in Groups and 
Clusters 



As discussed in M+10, we find strong evidence for a 
link between the cool core and the observed Ha filaments. 
Summarized briefiy, the main pieces of evidence in favor 
of this scenario are: 

1. The soft X-ray and Ha morphologies are correlated 
(Figure [l]) 

2. There is a correlation between the Ha luminosity 
and the X-ray cooling rate (Figure [5| 

3. Ha filaments extend all the way to the X-ray cool- 
ing radius, but never beyond (Figure |6|. 

4. The ICM coincident with the Ha filaments is cool- 
ing an order of magnitude faster than in the imme- 
diately surrounding regions (Figure [7|). 

We find that all of these trends, originally quoted in 
M+10, are seen in groups as well. This argues strongly 
in favor of a common origin and suggests that the cool- 
ing process may be detached from the global properties 
of the system. 

The link between the X-ray cooling properties and the 
optical emission suggests that the Ha filaments trace the 
cooling fiow. This is not to say that the Ha emission is 
due to the cooling ICM recombining at 10^ K, but rather 
that the morphology of the Ha filaments is reminiscent 
of magnetohydrodynamic simulations of gas cooling. As 
the ICM cools, it may be experiencing thermal insta- 
bility along magnetic field lines, causing it to collapse 
into thin, dense filaments resembling the observed Ha 
morphology (e.g., Sharma et al. 2010). Once the cool 
filaments have been established, they require an ioniza- 
tion source to produce the high Ha surface brightnesses 
which are observed. This heating may be due to a com- 
bination of photoionization by young stars, conduction 
with the hotter ICM and cosmic rays (see M+10 for a 
more detailed discussion). 

An alternative hypothesis is that of buoyant radio bub- 
bles (e.g., Reynolds et al. 2005, Revaz et al. 2008). In 
this scenario, bubbles are blown in the dense ICM by the 
radio-loud AGN. Due to their buoyancy these bubbles 
will be transported to larger radius, possibly entraining 
cool gas along the way. This model produces thin, radial 
filaments of cool gas reminiscient of the Ha filaments in 
many cluster cores. While radio- mode feedback is likely 
responsible for the optical emission in a few systems (e.g.. 
Hydra A, Abell 2052), it is hard to imagine how multiple 
filaments with similar surface brightnesses at a variety of 
orientations can be produced by this process. 

5. SUMMARY AND FUTURE PROSPECTS 

This study extends the work of McDonald et al. (2010) 
to include galaxy groups in order to determine the role 
of environment in the formation of Ha filaments in cool 
cores. In summary, we find: 

• The morphology and detection frequency of Ha fil- 
aments in groups is similar to those seen in clusters. 

• There is no obvious dependence between Ha emis- 
sion and the temperature of the system. There is 
a weak correlation between the Ha emission and 
the mass of the system, as measured at r25oo- This 
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weak correlation may be strengthened by the inclu- 
sion of high-mass, non-cooling clusters, which are 
underrepresented in our sample. 

• There is a weak correlation between the presence 
of Ha filaments and the global ICM gas fraction 
(51% confidence level) and a strong correlation 
with global entropy (89% confidence level), such 
that Ha filaments are more frequently detected in 
systems with low entropy and high gas fractions. 

• There is a correlation between the presence of 
Ha filaments and the properties of the cool core, 
namely the X-ray cooling rate (72% confidence 
level), the average temperature in the inner 100 
kpc (55% confidence level) and the entropy at a 
radius of 50 kpc (84%). The correlations between 
the Ha emission and the cooling rate and entropy in 
the cool core confirm previous results (Cavagnolo 
et al. 2008; McDonald et al. 2010) and favor the 
scenario where these thin filaments are the result 
of rapid, thermally-unstable cooling of the ICM 
(Sharma et al. 2010). 

• The most extended Ha filaments extend to the ra- 
dius at which the hot ICM is cooling in ^ 3Gyr, 
but not beyond. This suggests a natural timescale 
for the formation of filaments. 

• ICM cooling is enhanced by roughly an order of 
magitude in regions with Ha filaments, compared 
to surrounding regions at the same radius. 

• The cooling efficiency (cooling rate per unit gas 
mass) is higher in groups than in clusters, as is pre- 
dicted by cooling in the absence of feedback. This 
is a manifestation of the weak dependence of cool 
cores with mass (second bullet) coupled with the 
strong dependence of gas content on mass. This 



correlation would be further improved by the in- 
clusion of high-mass, non-cooling clusters, which 
are underrepresented in our our sample. 

• We confirm Sun's (2009b) finding that AGN in cool 
core groups have much lower radio luminosity than 
those in clusters. We show that this could explain 
the amount of feedback scaling with environment, 
since cool cores in groups can only survive if the 
amount of radio-mode feedback is small. We ob- 
serve three systems with radio luminosities higher 
than the amount needed to remove the cool core, all 
of which have very disrupted cool cores and devi- 
ate from the relationship found for groups/clusters 
with closed feedback loops. 

Our findings suggest that cool cores, and the optical 
line-emitting nebulae commonly associated with them, 
can form in groups and clusters over a large range of 
mass and temperature, with very little dependence on 
either. However, considerably larger samples will be 
needed to quantify whether this correlation is weak or 
nonexistent. We confirm the results of M+10, offering 
further support that the observed Ha filaments in the 
cores of groups/clusters are intimately linked to the X- 
ray cooling in most systems. Without high-resolution 
UV imaging and optical spectroscopy, we are unable to 
constrain the ionization mechanism of these filaments. 
We plan on addressing these issues in upcoming papers. 
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Table 1. - Properties of our sample of galaxy groups and clusters. X-ray temperature measurements (last coli 
- the superscript identifies the source: (1) Spatially averaged temperature (White et al. 1997), (2) Spatialb 

100 kpc (M+IO), and (3) Temperature at r25oo (Sun et al. 2009a). This sample covers a broad range in temperature, with values ranging 
from 1-10 keV. 
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Fig. 1. — X-ray and optical data for the 10 clusters in our sample. From left to right the panels are: 1) MMTF red continuum image, 
2) CXO X-ray image, 3) Unsharp masked CXO X-ray image, 4) MMTF continuum-subtracted Ha image. The horizontal scale bar in all 
panels represents 20 kpc. The X-ray and red continuum images are on the same scale, and the unsharp masked and Ha images are on the 
same zoomed-in scale. The square region in the X-ray panels represents the field of view for the zoomed-in panels. The grayscale in all 
images is arbitrarily chosen in order to enhance any morphological features. 
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Fig. 1. — Continued. 
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Fig. 2. — X-ray scaling relations for the 10 groups (blue) and 17 clusters (red) in our sample. The addition of groups to our sample allows 
us to probe the low temperature, mass and entropy regime. We differentiate between systems with Ho; emission (filled circles) and those 
without (open circles). Systems that are clearly disturbed either by mergers or AGN feedback are depicted by triangles and typically have 
under-estimated M25oo- The high-temperature, undisturbed outlier in panel (d) is Ophiuchus, which lies in the Galactic plane and, thus, 
suffers from heavy extinction. The high-entropy outlier in panels (b), (c) and (f) is Abell 2151, which requires the largest extrapolation 
from the last data point to r25oo, and thus has the least constrained properties. For comparison, we show relations derived from Sun 
et al. 2009a for several of the same systems (dashed lines). We show the functional form of these relations in the upper right corner. The 
scatter in these relations (dotted lines) represents our uncertainty in converting quantities from rsoo to r25oo- For panel (b), we plot the 
exact relation quoted by Sun et al. 2009a. 
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Fig. 3. — Frequency with which we observe Ho; emission in clusters and groups, as a function of the global X-ray properties. The black 
histograms show the frequency with which we detect any Ho; emission whatsoever, while the red histograms show the frequency with which 
we detect extended Ho; filaments. The horizontal lines represent the overall detection rate, while the vertical errorbars represent the 1-cr 
uncertainty in a given bin. The absolute number of systems in each bin are shown above the black histograms. We find that there is no 
measureable correlation between the global temperature and the presence of Ho; emission. There is a weak correlation between the presence 
of Ha and the global mass, and stronger correlations with both entropy and gas mass fraction, however these are at the < 1-cr level. The 
black and red numbers represent the F-test confidence with which we can reject the hypothesis of a steadily rising/falling distribution in 
favor of a flat distribution, for all Ho; emission and emission in fllaments, respectively. 
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Fig. 4. — Frequency with which we observe Ho; emission in clusters and groups, as a function of the X-ray properties in the group/cluster 
core. The black histograms show the frequency with which we detect any Ha emission whatsoever, while the red histograms show the 
frequency with which we detect extended Ho; filaments. The horizontal lines represent the overall detection rate, while the vertical errorbars 
represent the 1-a uncertainty in a given bin. The absolute number of systems in each bin are shown above the black histograms. We find 
that the correlation between Ha and 1.4 GHz flux to be slightly weaker than between the Ho; flux and the X-ray cooling rate. We also 
see anti-correlations between the presence of Ho; emission and the core temperature and entropy, as reported in M+10. The black and red 
numbers represent the F-test confldence with which we can reject the hypothesis of a steadily rising/falling distribution in favor of a flat 
distribution, for all Ho; emission and emission in fllaments, respectively. 
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Fig. 5. — Ho; luminosity, quoted as a star formation rate based on Kennicutt (1998), plotted against the average temperature in the inner 
100 kpc (kTioo), specific entropy at 30 kpc (K30), and the integrated mass deposition rate (dM/dt). Red points in this plot represent 
clusters from M+10, while the blue points refer to the groups sample in this paper. The SFRs in the lower panels have had the nuclear 
contribution removed (SFR*). The open star in the lower panels refers to Perseus A (Conselice et al. 2001, Sanders et al. 2004). The 
diagonal dashed line represents the limit where all of the cooling X-ray gas turns into stars, while the dotted line is the case where all of 
the cooling X-ray gas is made up of hydrogen which recombines only once (Fabian et al. 1984). 
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Fig. 6. — Left: Correlation of Kcooh the radius at which the coohng time of the ICM reaches 5Gyr; and Kna, the largest radius at which 
we detect Ha emission. When there was no detected Ho; emission an upper limit of the seeing FWHM has been enforced. The dashed line 
refers to the one-to-one case, while the point colors refer to the clusters from M+10 (red) and the groups from this paper (blue). There 
appears to be an upper limit on the radius of Ho; filaments, corresponding to the cooling radius. The black point represents a lower limit 
estimate on the cooling radius of Perseus A (Fabian et al. 2000, Conselice et al. 2001). Right: Similar to plot on the left, but using a 
cooling time of 3Gyr to define the cooling radius. Perseus, Abell 1795, and Sersic 159-03 lie almost exactly on the 1-to-l line in this case. 
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Fig. 7. — Upper panels: Distribution of temperature, electron density, entropy, and cooling time ratios in and out of the Ha filaments. 
The solid black and dashed grey lines bracket the extreme cases of the filament geometry: solid black lines are the thin-filament case, which 
is modeled with a two-temperature plasma, while the grey dashed lines are the case of single-temperature sheets of gas seen edge-on. The 
blue, dashed lines show the contribution to the total histogram from galaxy groups. The in-filament gas shares similar properties with the 
cool core, namely that the X-ray gas has a cooling time of ~ 10% that of the surrounding ICM. There appears to be no difference in the 
properties of the filaments in groups and clusters. Lower panels: Similar to the above panels, but now considering only the outer 70% of 
the filaments in radius, in order to remove the contribution from the central, circularly symmetric, region. 
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Fig. 8. — Upper left: X-ray cooling rate (dM/dt) as a function of M2500, the mass enclosed within r25oo- Red and blue points refer to the 
cluster and group samples, respectively. Green points are six new groups taken from the Sun et al. (2009a) sample in order to better sample 
the X-ray properties of the groups. Upper right: X-ray cooling rate per unit total mass as a function of M25oo- We note an apparent upper 
limit on this scaled cooling rate at ~ 2 x 10~^ Gyr"-*^. Lower left: ICM cooling efficiency (cool rate per unit gas mass) as a function of 
^2500- The dashed lines reflect the expected slopes for Bremsstrahlung (-0.3) and line cooling (-0.96). The zeropoints for these lines were 
adjusted to best flt the data. The fact that the data matches the general predicted trend suggests that the effect of feedback must either 
be negligible or scale with total mass, otherwise the trend would be altered. Lower right: Ha-determined star formation rate per unit gas 
mass versus M25oo- Similar trends are noted here since Lna is an excellent proxy for dM/dt (M+10). 
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Fig. 9. — 1.4 GHz radio luminosity versus integrated soft X-ray luminosity within the cooling radius. Color refers to the system mass, 
with blue, green and red referring to systems with M2500 ~ lO"*^^ M©, IQ-*^^ M©, and IQ-*^^-^ M©, respectively. Stars, closed circles and 
open circles are systems with filamentary, nuclear or no Ho; emission. Grey points are groups and clusters from Sun (2009b) and the red 
ellipses describe the two loci identified by Sun. The grey arrows show the proposed evolution of systems from cool, feedback-regulated cores 
(right ellipse) to X-ray faint corona (left ellipse). If the radio luminosity exceeds a certain threshold, depicted by horizontal dotted lines, 
the amount of feedback exceeds the required pdV work needed to disrupt the cool core. This simple picture seems to have some merit, as 
the three clusters in our sample with extremely disrupted cores (Hydra A, Abell 2052, Abell 4059) all have radio luminosities exceeding 
the allowable threshold for their mass and are significantly offset to the left of the cool core locus. This plot provides an explanation for 
the fact that groups with cool cores tend to have low-luminosity AGN when compared to their high-mass counterparts. 



